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FORMULAS FOR THE LAPLACE TRANSFORM OF 
STOPPING TIMES BASED ON DRAWDOWNS AND 

DRAWUPS 

By Hongzhong Zhang^ and Olympia HADJILIADIS*'t 
Graduate Center^ and Brooklyn College^, C.U.N.Y. 

In this work we study drawdowns and drawups of general diffu- 
sion processes. The drawdown process is defined as the current drop 
of the process from its running maximum, while the drawup process is 
defined as the current increase over its running minimum. The draw- 
down and the drawup are the first hitting times of the drawdown and 
the drawup processes respectively. In particular, we derive a closed- 
form formula for the Laplace transform of the probability density of 
the drawdown of a units when it precedes the drawup of b units. We 
then separately consider the special case of drifted Brownian motion, 
for which we derive a closed form formula for the above-mentioned 
density by inverting the Laplace transform. Finally, we apply the re- 
sults to a problem of interest in financial risk-management and to the 
problem of transient signal detection and identification of two-sided 
changes in the drift of general diffusion processes. 

1. Introduction. In this paper we derive the Laplace transform of the 
probability density of the drawdown of a units when it precedes the drawup 
of b units for a general diffusion process. The drawdown process is defined 
as the current drop of the process from its running maximum, while the 
drawup process is defined as the current increase over its running mini- 
mum. The drawdown and the drawup are then the first hitting times of the 
drawdown and the drawup processes respectively. The derivation is first ac- 
complished in the case that a = b, by drawing the connection of the relevant 
event to the range process. We then consider the case a 7^ b and derive the 
Laplace transform through path decomposition. A key element in the above 
derivation is a function related to the first hitting times of the underlying 
diffusion process. We then consider the special case of drifted Brownian mo- 
tion with drift ji and volatility a. In this case we are able to invert the 
Laplace transform and derive the density of the drawdown of a units when 
it precedes the drawup of b units, which we denote by p^(t;a,b). Finally, 
we discuss the applications of the results to a problem of interest in finan- 
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cial risk-management and to the problem of transient signal detection and 
identification of two-sided changes in the drift of general diffusion processes. 

Our results extend the work of Taylor [23] and Lehoczky [13]. Taylor 
[23] derives the joint Laplace transform of the drawdown and the maximum 
stopped at the drawdown in a drifted Brownian motion model. Lehoczky 
[13] extends the above result in the case of a diffusion process. However, 
none of these results relate the drawdown to the drawup. In a recent paper 
of Salminen & Vallois [21], the joint distribution of the maximum drawdown 
and the maximum drawup processes is studied in a drifted Brownian motion 
model; yet it is not possible to extract information on the joint distribution 
of the drawdown and the drawup from this paper. In Hadjiliadis & Vecer 
[10], a closed- form formula is derived for the probability that the drawdown 
precedes the drawup in a drifted Brownian motion model. This result is later 
extended to diffusion processes in Pospisil, Vecer Sz Hadjiliadis [19]. In Zhang 
& Hadjiliadis [27], the authors obtain the probability that the drawup of a 
units precedes the drawdown of equal units in a drifted Brownian motion 
model in a finite time-horizon. However, the approach used there only ap- 
plies to a drifted Brownian motion model, cannot be extended to a general 
diffusion process, and does not work in the general case a ^ b. In this paper, 
we supersede these limitations through the Laplace transform method. 

Drawdown processes have been extensively used in the financial risk- 
management literature. Grossman &; Zhou [7], Cvitanic & Karatzas [5], 
Chekhlov, Uryasev & Zabarankin [4] studied portfolio optimization under 
constraints on the drawdown process. Magdon-Ismail et. al. [14] determined 
the distribution of the maximum drawdown process of Brownian motion, 
based on which they described another time-adjusted measure of perfor- 
mance known as the the Calmar ratio (see [15]). Meilijson [16] proved that 
the drawdown can be viewed as the optimal exercise time of a certain type 
of look-back American put option. Other works which describe drawdown 
processes as dynamic measures of risk include Vecer [24, 25], Pospisil & Ve- 
cer [18], Pospisil, Vecer & Xu [20], Zhang & Hadjiliadis [27]. An overview of 
the existing techniques for analysis of market crashes as well as a collection 
of empirical studies of the drawdown process and the maximum drawdown 
process please refer to Sornette [22]. 

Drawdown processes do not only provide dynamic measures of risk, but 
can also be viewed as measures of "relative regret". Similarly drawup pro- 
cesses can be viewed as measures of "relative satisfaction" . Thus the a draw- 
down or a drawup of a certain number of units may signal the time in which 
an investor may choose to change his/her investment position depending 
on his/her perception of future moves of the market and his/her risk aver- 
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sion. Using the results in our paper we are able to calculate the probability 
that a relative drawdown of (100 x a)% occurs before a relative drawup of 
(100 x f3)% in a finite time-horizon. On the other hand, a digital option 
on the event that the relative drawdown occurs before the relative drawup 
could also be seen as a means of protection. Our paper provides a closed- 
form formula for the risk-neutral price of this digital option at time both 
in the case of an infinite maturity and in the case of a finite maturity. 

Drawdown and drawup processes have also been used in the problem of 
quickest detection of abrupt changes in a stochastic process. More specifi- 
cally, consider the situation in which a diffusion process is sequentially ob- 
served. At some unknown point in time, possibly as a result of the onset of 
a signal, the dynamics of the process change abruptly in one of two possi- 
ble opposite directions in the drift. Drawdowns and drawups then provide 
a detection mechanism of the change point for each of the possible changes. 
They are known as CUSUM stopping times in the statistics literature, and 
their properties have been extensively studied (see Barnard [1], Dobben [6], 
Bissell [2], Woodall [26], Hadjiliadis & Moustakides [9], Khan [12], and Poor 
& Hadjiliadis [17]). A challenging problem in engineering is the detection 
and identification of such signals when they are only present for a finite pe- 
riod of time. These signals are known as transient signals. Using the results 
in this paper, it is possible to derive closed-form formulas for the probabil- 
ity of misidentification of the direction of the change in the drift when the 
signal has exponential life. Moreover, using our results for drifted Brownian 
motion, we derive this probability when the transient signal is present for a 
finite period of time T. 

The paper is structured in the following way: definitions and a fundamen- 
tal lemma are introduced in Section 2. In Section 3, we derive the Laplace 
transform of the drawdown of a units when it precedes the drawup of b units, 
in the cases a = b (Theorem 3.1), a > b (Theorem 3.2) and a < b (Theorem 
3.3). A special case of a drifted Brownian motion model is discussed in Sec- 
tion 4, where we also derive the closed-form density p^\t; a, b) by analytical 
inversion of the Laplace transform. We then present applications of our re- 
sults in a problem of risk-management and the problem of transient signal 
detection and identification of two-sided alternatives in Section 5. Finally, 
we conclude with some closing remarks in Section 6. 

2. Drawdown and Drawup processes. We begin with a mathemat- 
ical definition of a drawdown and a drawup in a diffusion model and present 
the fundamental lemma. 

Consider an interval / = (l,r), where — oo < / < r < oo. Let (0,,^, P) 
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be a probability space, {Bt,t > 0} a Brownian motion, and {Xt,t > 0} the 
unique strong solution of the following stochastic differential equation: 

(2.1) dX t = fx{X t )dt + a(X t )dB t , Xq = x E I, 

where Xt E / for all t > 0, and cr(u) > for all u S I. We will use P x (-) to 
denote P(-\Xq = x). 

The drawdown and drawup processes are defined respectively as 

(2.2) DD t = sup X s - X t , 

0<s<t 

(2.3) DU t = X t - inf X s . 

0<s<t 

The drawdown of a units and the drawup of b units are then defined respec- 
tively as 

(2.4) T D (a) = inf{f > 0\DD t =a}, a> 0, 

(2.5) Tu(b) = inf{i > 0|DJ7 t = 6}, b > 0, 

where, by convention, we assume that inf = oo. 

In the following section, we derive the main results in this paper. We need 
the following fundamental lemma to finish the proofs. 

Lemma 1. Let us denote by t u , u £ I , the first hitting time of the process 
{X t ;t> 0} to u. That is, 

t u = mf{t > 0\X t = «}, 

For y < x < z and A > 0, define 

(2.6) E x [e~ Xr y • l {Ts<T , } ] := i X {y,z;x,\). 
Then 

(9 7\ P x (~, \\- 9{x;\)h{z;\) - g(z; X)h{x; X) 

[ ] ~ d(y, X)h(z; A) - g(z; X)h(y; A) 

with g(-;\) and h(-;X) being any two independent solutions of the ordinary 
differential equation 

(2.8) f'W^^wf V 

Proof. See [13], page 603. □ 
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3. Mathematical Results. In this section we derive formulas for the 
Laplace transform of the probability density of the drawdown of a units when 
it precedes the drawup of b units, for any a, b > satisfying x ± a, x ± b G 
I = [I, r]. We have 



(3.1) E x 



-XT D {a) -, 
e ■ 1 {T D (a)<T u (b)} 



e- xt P x (T D (a) € dt,Tu(b)>t). 



In the sequel we denote the Laplace transform (3.1) by (A; a, b). 

3.1. The case of a = b > 0. We determine L x (X;a, a) for a > in this 
paragraph. 

Theorem 3.1. For a > and A > 0, we have 

(3.2) L^(\;a,a)= \ —£ X (u,u + a;x,\)du. 

Jx-a oa 

Proof. First, it is easily seen that for t > and a > 0, 

(3.3) {T D {a) G dt,Tu(a) > t} = {p{a) e dt,x - a < X t < x}, 
where p(a) is the first hitting time of the range process 

(3.4) p(a) = inf{t > 0\R t = a}, 
with 

(3.5) R t := supXs-infX, = DU t + DD t . 

s<t s <t 



So it suffices to determine 



E , 



e- xt ■ 1 



' L {T D (a)&dt,T u {a)>t,X t =u} 



for all x — a < u < x. For this purpose observe that for any t > 0, a > and 

x — a < u < x, 

{Td(o) £ dt,Tjj(a) > t,X t = u} = {r u G d£,supX s = u + a}, 

s<t 
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which suggests that for any A > 0, 



E 



-Alb (a) . i 



{T D (a)<Tu(a),X TD{a) =u} 



-XI 



-XI 



■ E_ r . 



■ Ex 



1 {T D (a)Gdt,T (7 (a)>i,Js: i =n} 



l{i~ u e<it,sup s<t X a =u+a} 



Ft 

e~ xt ■ ^E x 
oa 

Oa 



^-{Tv.edt, sup s<t X s <u+a} 



{T u =dt,T u + a >t} 



Oa 



p -Xt u -i 

- L {T-u<T„+a} 



From Lemma 1 it follows that 



(3.6) E x 



-XT D (a) . 1 



{T D (a)<T c/ (a),X T£l(a) =u} 



d_ 

da 



£ x (u, u + a;x, A). 



Then integration of the above identity over the interval (x — a, a) in u yields 
(3.2) and completes the proof of the theorem. □ 

3.2. The case of a > b > 0. We determine L x (X;a,b) for a > b > in 
this paragraph. To prove the main result we need following proposition. 

Proposition 1. For b > 0, c < u satisfying c,u + b G I, and A > 0, 
define 



H x (X;b,c) 



(3.7) 
Then 

(3.8) H*(\;b,c) 



En, 



— Xt c -I 
e • l {Tc<Tu(b)} 



exp 



« d 
dw 



l x (v,v + b;w,X)dv 



Proof. We follow the idea of [8], [13], , and partition the interval [c, u] 
into k subintervals {[vi, 1 < i < k} with u = vo > v\ > . . . > vj~ = c. 
Let Afc = maxi<j</ c (wj_i — i>j) and assume — > as k — > oo. As a discrete 
approximation to H x (X;b,c) defined by (3.7), compute 



E u e - A EtiK- r ^-i) . i 

k 



i=i 



e" Ar ^ • 1 



{after r Vi _ 1 ,Xt hits «i before increasing to Vi+b,l<i<k} 
{r Vi <T v , +b } ) 
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where the last equality follows from the strong Markov property. 

It will be shown that as k — > oo and A& — > 0, the limit of above expression 
exists and does not depend on the particular sequence of partition chosen 
and hence is equal to H x (A; b, c). 
By Lemma 1, 



i=i 

Taking log gives us 
k 



{r Vi <T vi+b } 



~[£ x (vi,Vi + b;vi-i,\). 



i=l 



J^log t (vi,Vi + b; Vi-i,X) 

i=l 
k 

^ log £ x (vi,Vi + b;vi, A) + £ x (v,i,Vi + b;i)j_i,A) - £ x (vi,Vi + A) 



A' 



5^ log 1 +£ x (vi,Vi + b;Vi-i,X) -£ x (vi,Vi + b;vi,\) 




i=i 

k 

E 

i=l 

U Q 

dw 



dw 



£ x {vi,vi + b-w,\) ■ (vi-! - Vi) + 0(A fc ) 



+ &;tu, A)cfo, as A fc -> H 



from which we obtain 

H x (X;b,c) = exp 



ii ^ 
dw 



£ x (v,v + b;w,X)dv 



This completes the proof of Proposition 1. 

Now let us state and prove the main result in this paragraph. 
Theorem 3.2. For a > b > and X > 0, we have 

(3.9) L x (X;a,b) = / — £ x (u, u + b; x, X) ■ H x (X; b, u - a + b)du. 

J x—b OO 

Proof. Any path in the event {To (a) < Tjj(b)} has the decomposition 

1. {X t ;0<t<T D (b)}; 

2. {X t+TD{b) ;0< t < T D (a) - T D (b)}. 



□ 
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Conditioning on Xx D {b) = u -> the process in 2. starts at u, and decreases to 
u — a + b before it incurs the drawup of b units occurs. This gives rise to the 
representation 

(3.10) T D (a) = T D (b) + r u . a+b o 9(T D (b)). 

Therefore, for x — b < u < x, 

-\T D (a) i 
e i {T £) (a)<T c/ (6),X T£)(6) = U } 

~ 6 1 {T D (a)<T u (b),X TDW =u} 
— c -AT D (fe) -, 



(3.11) 



before T D {b) 
v „-Ar u _ a+ f,o6»(T (b)) 

x e a+b . l {Tti _ a+b0e(TD(fe))<Tc/(fe)oe(Tl3(6 ) )} 

after T D (6) 



To get the expectation of the above expression under E x [-] } we first compute 
its conditional expectation under £' u [-|^ r r£) ( 6 )]. By the strong Markov prop- 
erty, the factor "before Td(£>)" is deterministic under E u [-\J r TD ^ : and the 
factor "after Tb (6)" has conditional expectation 



E a 



1 K_ a+6 °0(T D (&))<T t7 (&)o0(T J3 (&))} I >T D (6) 
{T u -a+b<Tv(b)} 



e —\r u - a+b _ 2 



which, by Proposition 1, is equal to H^(X; b,u — a + b). Taking the expec- 
tation of (3.11) under -E x [ - ]j and using (3.6), we obtain 



Er. 



E, 



-XT D (a) i 
e L {T D {a)<Tu{b),X TD(h) =u} 

e -AT D (6) . l {TD{b)<Tu ^ XTD(h _ u} ] .H*(\;b,u-a+ b) 



(3.12) = ^e x (u,u + b;x,X) -H* (X;b,u - a + b). 

The integration of the above identity over the interval (x — b, x) in u yields 
(3.9) and completes the proof of Theorem 3.2. □ 

3.3. The case of b > a > 0. We determine L^(X; a, b) for b > a > in 
this paragraph. To prove the main result we need the following proposition. 
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Proposition 2. For any a > and x £ I satisfying x — a 6 I , and 
A > 0, define 



(3.13) 

Then 

4? (A; a) 



J* (A; a) 



-AT D (a) 



yjX/" \\ — f U jt-\ £ x (v,v — a:w,X)dv , 

i (u — a, u; w, A) ■ e J * 9w\ w=v ; q^. 



Proof. See [13], page 602. 

Now let us state and prove the main result in this paragraph. 
Theorem 3.3. For b > a > and A > 0, we have 



x+a pi 

- ° -2x-X 



□ 



dv—t x ~* (2x - v, 2x - v + a: x, A) 
da 



L*(A;a,6) = J*(A;a)- 
(3.14) xH%l?(\; a,2x-v-b + a)- j£ 6 _ a (A; a) 

Proof. First, it is easily seen that for b > a > 0, 



L*(\;a,b) = J£(\;a)-E : 



X, 



e -XT D (a) . 1 



{sup s < Tl)(a) DC/ s >fe} 



Therefore, to prove (3.14), it suffices to show that 



Em 



e -\T D (a) . x 



( su P S <T D {a) DU s >b) 



:+a Q 

dv—£ 2x ~ x (2x - v, 2x - v + a; x, 
da 



X 

(3.15) xF^lf (A; a, 2x - v - b + a) • J,f +6 _ (A; a 

To do so we observe that for b > a, 



(3.16)^ 



-ATo(a) . ]_ 



-ATo(a) . 1 



{Tu(b)<T D (a)} 



And we are going to calculate the right hand side of (3.16) in the sequel. 

Consider the path decomposition for any path in the event {Tjj(b) < 
T D (a)}. We have 

1. {X t ;0<t<Tu(b)}; 

2. {X t+Tu{b) ;0 < t < T D (a)-Tu(b)}. 
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Intuitively, before time Tjj(b), the process experiences no drawdown of a 
units and the first drawup of b units occurs at Tjj(b), when the process also 
reaches a new maximum; thereafter, the process has a drawdown of a units 
at time To (a). Thus for any path in the event {Tjj(b) < To (a)} we have 

(3.17) T D (a) = T ir (b) + T D (a) o 9{Tu(b)). 

Therefore, for b > a and x < v < x + a, 



(3.18) 



E x 



e -XT D (a) . h 



{sup t <T D ( Q) 

e -AT c/ (6)+T £) (a)o0(T u (6)) . j 



DU t >b,X Tu(a) =v} 

{ su Pt<T D (a) DU t >b,X Tu(a) =v] 



E x 



E x 



E x 



p-ATtr(6) . i _ _ ^ . ^ v , x ATq (a)o9(Tu (b)) 



after Tj/(b) 
-AT D (a) 



before 2^(6) 

• i {T c/ (6)<T (a),X T[/(a) = t ;}^+6- 

" AT t/( 6 ) 1 7 X 

1 {r C /(fe)<T c (a),X T[/((l) =^} • ^+6-0 (A «J 

• i {T [/ (6)<T D (a),X Tc/((l) =^}J • J„ +6 _ a (A,aj, 



where the third equality follows from the strong Markov property. The ex- 
pectation in the last line can be determined as follows. Note that for the 
process {Y t = 2x — X t , t > 0}, 

dY t = -p{2x - Y t )dt + a(2x - Y t )dB' t , Y = x 

with B t = —B t , the vector of random variables (T^ (6), T^f (a), 2x — Yrpy^-^j 1 

has the same law as the vector of random variables \ TD(a),Tu(b),X Tu ( a 
for {Xt;t > 0} under ^[-J. So we know from (3.12) that 



(3.19) 



E x 

E x 

d_ 

Oa 



e -XT v (b) -, 

e • 1 {T u (b)<T D (a),X Tu(a) =v} 



l {T Y {b)<T^{,a)X rA 



Tlia)= 2X -^ 



iovX t 
for Y t 



£ 2x ~ x (2x - v, 2x - v + a; x, A) • H'glf (A; a, 2x - v - b + a). 



2x-X, 



The integration of (3.18) over the interval (x,x + a) in v yields (3.15) and 
completes the proof of Theorem 3.3. □ 

We now proceed to treat separately the case of a drifted Brownian motion 
model. 



1 T D (6) and T v (a) are the drawdown and drawup of the process {Yt; t > 0} respectively. 
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4. A study of Brownian motion. In this section, we apply the results 
in Theorem 3.1, Theorem 3.2 and Theorem 3.3 to a drifted Brownian motion 
model and calculate the probability density of the drawdown of a units when 
it precedes the drawup of b units. 

First it is easily seen that / = (— oo, oo) in a drifted Brownian motion 
model. The function £ x (y, z; x, A) for Xt = x + /it + crWt can be found in [3] 
page 295: 

(4-1) ^(y,z;x,A) = 

smh[(z - y)5 M)CT (A)J 

where S flj(T (X) = y / (2\/a' 2 ) + (// 2 /<t 4 ). Thus the Laplace transforms in The- 
orem 3.1, Theorem 3.2 and Theorem 3.3 can be calculated explicitly as: 
1. a = b > 0: 

T x n \ _ S^(X) Je-^[^ )CT (A)coth[a^, CT (A)] + ^] 
L [A; a) — 



(4.2) 

2. a > b > 0: 




sinh[a5 , Mi(T (A)] 



sinh [aS M)CT (A)] 



(4.3) Lq (A; a, b) = Lq (A; b) • exp [T^A; 6) (a - 6)] , 
where 

(4.4) r M , CT (A; b) = -^- S^ a {\) coth[6^ )CT (A)]; 
3. b > a > 0: 

(4.5) L^(A;a,6)= [l-L - x (A;a)-e T --( A ^- a )] • J x (A;a), 
where 



(4.6) 



S^ a {\) cosh[aS Mi(T (A)] - (/i/a 2 ) sinh[a5 Ali<J (A)] 
S' Mi(7 (A;a)e~^ 1 



sinh[aS Mj(J (A;a)] T_ Mj(T (A;a) 

One can easily obtain several known results from (4.2), (4.3) and (4.5). 
First, by letting A — > + , the formulas coincide with the probability results 
in Hadjiliadis & Vecer [10]. Second, by letting b — > oo in (4.5), one obtains 
the Laplace transform of Tjj(a), J^(X; a). 

Moreover, we can invert (4.3) analytically to obtain the density P(Td(o) G 
dt, Tjj{b) > t) for any a > b > 0. In fact we have 
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Theorem 4.1. Define p^(t; a, b)dt = P{T D (a) G dt,T v {b) > t) for 
a > b > 0, then 



pM(t;a,b) 



2 M 2 t M(a-b) 

— e 2ct 2 a 5 



E 

m,n=0 



(m + n + 1)! 
(m + l)\m\n\ 



2(a - b) 
ay/I 



2fi 2 »(a-6) ~ ( m + n +l)l 



-e ^ 



E 

m.n=0 



(m + l)\m\n\ 



2/i(a - b) 



0" 



(4.7) 



-l) m e-^G^(t) 



where 5 n is the Kronecker delta and 

2k 



fc=0 
m+1 



(7 



(m+l-2fc) 



(2m + 2n + 1)6 + a 



cry 7 * 



fc=0 
m+] 

^ 3) (*) = E 



fc=0 



fc r- 



1 + (-!)' 



. m + n + 2 



j . /i(2mb+a) 



( (m+l-fe) 

2to6 + a + /it 



n + 1 
2m6 + a 



(m+l-fe) 



2(m + n + 1)6 + a 



a 



Vt 



a 



Vt 



ay/i 



with <fi and <3? being the standard normal probability density and cumulative 
distribution respectively. 4>^ is the k-th derivative of <f>. 

Proof. We start by rewriting (4.3) in a more tractable way, 
(aq\[°j i^5' Ati(J (A)sinh[(-n)5 / , i(T (A)] 

m j_ b due^ smh2[5 c^ (A)] exp[T At , CT( A;6)(a-6 ) ]. 

By using the first formula on page 643 of [3], in their notation, we obtain 
the inverse Laplace transform of the integrand in (4.8) 

_ e ^ [es a 2 t (l,2,b,u,a - 6) - es CT 2 t (l, 2, 6, -u,a - 6)] . 
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After some simple manipulation, the above expression becomes 



(4.9) 



2e 2a 



uTt | ji{u-\-b — a) 



aVt 3 

(2m + 2n + 1)6 + a + u~ 



E 



(m + n + 1)! 



2(a - b) 



ay/i 



a-s/i 



_ (m + l)!m!n! 

m,n=0 v ' 

(2m + 2n + l)b + a-u 



a 



Formula (4.7) follows from integration of (4.9) over (—6, 0) in u. 



□ 



One can let a = b in (4.9) to get a similar joint probability density as that 
in Proposition 1 of [27]. 

Similarly, we can invert (4.5) analytically to get the joint density 



P(T D (a) G dt, sup DU S G a + dz) 

8<t 

In particular we have 



d 

—p^(t;a,a + z)dtdz Va,z>0. 
oz 



Theorem 4.2. For a Brownian motion with constant drift fx and con- 
stant volatility a, any a, z > 0, we have 

fj, 2 t | n(z-a) 



^pM(t;a,a + z) 



4e 2a 



E 



(m + n + 2)! / 2z 



\ m 



: 



aV¥ J^ =Q {m + 2)\m\n\ \aV~tJ 

2F$> n (t, z) - e^F$ n (t, z ) - 5 n e^F^(t, 
4^3 „ (z _ a) ~ ( m + n + 2)i /2fxz\' 
~a* e ° ^ n (m + 2)!m!ra! \ 

m,n=0 v ' / x 



G 



(/<) 



( M ) _ (-lrG^ i(t,ar) - e^G^ +n (t,z) 



(4.10) 

where 

-^m,n(^) ^) 



+( -l)- e 5 G ^ +1 (t,^ 



L^J r 

: E 

fc=0 - 



2A- 



(m+2-2fc) 



(2m + 2n + 3)a + z 



k=0 
m+2 

*1 3) (*,*) = £ 



a 
a 



(_1)* + m + n + 3 



n + 1 



cr\/i 

,( m+ 2-fc) 



(2m + 2n + 4) a + z 



a 



(m+2-k) 



k=0 

M[2(m+l)a + z] 



(2m + 2)a + z 



2(m +l)a + z + fit 
ay/t 
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Proof. We start from the equality 
-- J X (X; a) - Lq x {\; a ) e T-,A^){b-a) jX ^. a) 

-J x (A; a) - L^ x (X;a)J x (X;a) + S^;a)Lo X (X;a) ^^^^ 

sinh[aS' Atj0 -(A; a)]e^ •'O 
[S„ ff (A; a)] 2 sinh[(-n) < S M , (T (A; a)] 



/U ro—a 
du / tfo^ 
-a 7o 



sinh [aS^alX; a)\e <?' 1 



By using the first formula on page 643 of [3], the integrand in the last line 
has inverse Laplace transform 

—e *Z* T 2 [es a 2 t (2,3,a,u,z) - es a 2 t (2,3,a,-u,z)]. 
After some simple manipulation, the above expression becomes 

4e~^~ ^ (m + n + 2)\ f 2z 



^(m+3) (2m + 2n + 3)a + g - u _ (m+3) 
cryt 



(2m + 2n + 3)a + z + it 



x 



The integration of (4.11) over (— a, 0) in it yields (4.10) and completes the 
proof. □ 

5. Application. 

5.1. Relative drawdowns and relative drawups of stock prices. Consider 
the case of a stock with geometric Brownian motion dynamics under a prob- 
ability measure P: 

(5.1) dS t = fiStdt + aStdWt, Sq = 1. 

Using Theorem 4.1 and Theorem 4.2, we are in the position to address 
the following question: 

What is the probability that this stock would drop by (100 X ct)% 
from its historical high before it incurs a rise of (100 X /3)% from 
its historical low in a pre-specified plan horizon T? 

First observe that 

(5.2) dlogS t = udt + adWt, log So = 0, 
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where v = /j — \a 2 represents the logarithm of the return of the stock. 

Now define the running maximum and the running minimum of the stock 
process {St} 

M t = sup S s , N t = inf S s . 

s <t s <* 

We also let Ud{oi) be the first time the stock drops by (100 x a)% from its 
historical high and Ur((3) the first time that the stock rises by an amount 
equal to (100 x /?)% from its historical low. That is, 

(5.3) U D (a) = M{t>0\S t = (l-a)M t }, 

(5.4) U R ((3) = M{t>0\S t = (l + (3)N t }. 

Thus, it is possible to calculate the exact expression for the probability 
that a percentage relative drop of (100 x a)% precedes a relative rise of 
(100 x 0)% by noticing that 



(5.5) 



U D (a)= T D (-log(l-Q)) 
Ur(P)= T v Qog(l + P)) 



And this probability can be calculated explicitly as 

P(U D (a)AT<U R (p)AT) = /%<">(*; - log(l - a), log(l + (3))dt. 

Jo 

Moreover, a digital option on the event that the relative drawdown pre- 
cedes the relative drawup can also be perceived as a means of protection 
against adverse movements in the market. In particular, the discounted pay- 
off of this digital option can be written as 

(5.6) PO(a,f3) = e~ rt • l{u D ( a )edt,u R (p)>t} • 1 {t<T}, 

where r > is the risk- free interest rate and T is the maturity of the option. 

Under a risk-neutral measure Q, the stock price and its logarithm have 
the following dynamics respectively, 

(5.7) dS t = rS t dt + aS t dW t , S = 1, 

(5.8) dlogS t = v'dt + adWt, log So = 0, 

where v = r — \(J 2 ■ 

Using (5.5) and our results we are able to derive the risk-neutral price at 
time of this digital option: 
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/ e~ rt 


[j£ 


/o 





In the case of a perpetual option (see [11]), the risk- neutral price of the 
digital option is already given by the Laplace transform (4.2), (4.3) and 
(4.5). In particular, 

Q[PO(a,P)] = 4 og5 (r;-log(l-a),log(l+/3)). 

In the case of a finite life option maturing at time T < oo, we can apply the 
densities (4.7) and (4.10) to calculate the risk-neutral price. 

1. (l-a)(l+/3) < 1: 

(5.9) Q[PO(a,0)] = [ T e- rt p^\t;-log(l-a),log(l + P))dt; 

Jo 

2. 5 = (l-a)(l + /3) > 1: 

(5.10) Q[PO(a, p)] - Q[PO(a, a/(l - a))] 
rlog<5 q , 

—p {u >(t; -log(l - a),-log(l - a) + z)dz dt. 

5.2. The problem of transient signal detection and identification of two 
sided changes. In this example, we present the problem of transient signal 
detection and identification of two-sided changes in the drift of a general 
diffusion process. More specifically, we give the formulas for the probabilities 
of misidentifying the direction of the signal both in the case of exponential 
life transient signals and in the case of deterministic life transient signals. 
In particular, let {Xt,t > 0} be a diffusion process with the initial value 
Xq = x and the following dynamics up to a deterministic time r: 

(5.11) dX t = a{X t )dW u t < r. 

For t > t, the process evolves according to one of the following stochastic 
differential equations: 

(5.12) dX t = n(X t )dt + a(X t )dW t t > r, 

(5.13) dX t = -^{X t )dt + a{X t )dW t t > r, 

with initial condition y = X T . 

The time of the regime change, r, is deterministic but unknown. We 
observe the process {X t ,t > 0} sequentially and our goal is to detect the 
time of onset of the signal, as well as possibly identity its direction, before 
the signal disappears. 
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Let us denote by P£' + and P%~ the probability measures generated on 
the space of continuous functions C[0, oo) by the process {X t , t > 0}, if the 
regime changes at time r from (5.11) to (5.12) and from (5.11) to (5.13), 
respectively. 

In this context suppose that the drawdown of a units, To(a), can be used 
as a means of detecting the change in the dynamics of {Xf,t > 0} from 
(5.11) to (5.13). Then the drawup of b units, Tjj(b) may be used as a means 
of detecting the change in the dynamics of {Xt, t > 0} from (5.11) to (5.12). 
For example, in the case that //(•) = fi > 0, it is easy to see that the drawup 
could be used as a means of detecting a change from (5.11) to (5.12), while 
the drawdown would provide a means of detecting a change from (5.11) to 
(5.13). In particular, T(a,b) = mm{Tu(a),Tu(b)}, also known as the two- 
sided CUSUM (see Khan [12]), has been extensively used as a means of 
detecting two-sided alternatives in the drift (see Barnard [1], Dobben [6], 
Bissell [2], Woodall [26], and Poor & Hadjiliadis [17]). In the special case in 
which /z(-) = fx > and there is no reason to believe that a change from 
(5.11) to (5.13) is more or less likely than a change from (5.11) to (5.12), 
it is natural to use thresholds a = b. However, in the general case different 
thresholds a and b could be used. 

In many applications in engineering, the life of the signal after its onset is 
often limited. The lifetime of the signal is also random and may not depend 
on the dynamics of the underlying process. In the case of exponential life 
transient signals, the signals are present (after r) for a period of time £, 
where C is an independent exponentially distributed random variable with 
parameter A > 0. Theorems 3.1, 3.2 and 3.3 can be used to compute the 
probability of sequential misidentification of the signal in the case that the 
onset of the signal occurs at time 0. More specifically, 



P%' + {T D (a) AC < Tu(b) AC) = / Pq {Td(cl) A i < Tjj(b) A t) • Xe dt 



where follows (5.12) with r = 0, expresses the probability that an alarm 
indicating that the regime switched to (5.13) will occur before C and an alarm 
indicating that the regime switched to (5.13) while in fact (5.12) is the true 
regime. Thus (5.14) can be seen as the probability of a misidentification. 
Moreover, in the case that the density of the random variable X T admits a 




(5.14) 



K ' (A; a, b), 
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closed-form representation, we can also compute 

(5.15) J P T y + (T D {a) o 9(t) A C < T v (b) o 9{t) A Ofx T (v\x)dy 

= J Ly°' + (X,a,b)fx T (y\ x )dy, 

which can be interpreted as the aggregate probability (or unconditional 
probability) of a misidentification for any given change-point r. 

In the case of deterministic life transient signals, the signals are present 
(after r) for a finite period of time T. Using Theorem 4.1 we are still able 
to compute the probability of misidentification for Brownian motion (cr(-) = 
a > 0,/i(-) = fi). More specifically, 

(5.16) P^ + (T D (a) o9(t) AT < Tu(a) o 9{t) A T) = f p {ll) (t;a,a)dt, 

Jo 

expresses the probability of misidentification for any given change-point r. 

6. Conclusion. In this paper we derive a closed-form expression for the 
Laplace transform of the drawdown of a units when it precedes the drawup 
of b units for a general diffusion process. We then derive the probability den- 
sity of a drawdown when it precedes a drawup in the special case of a drifted 
Brownian motion model by inverting the Laplace transform. Although sev- 
eral authors in the literature have studied drawdowns and drawups ([23], 
[13], [10], [19], [21], [27]), this paper summarizes the probabilistic properties 
of a drawdown on the event that it precedes a drawup for a general diffusion 
process. These results are of practical interest in two main areas: financial 
risk-management and transient signal detection and identification. 
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